Mesoscopic conductance fluctuations in InAs nanowire-based SNS junctions 
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We report a systematic experimental study of mesoscopic conductance fluctuations in supercon- 
ductor/normal/superconductor (SNS) devices Nb/InAs-nanowire/Nb. These fluctuations far exceed 
their value in the normal state and strongly depend on temperature even in the low-temperature 
regime. This dependence is attributed to high sensitivity of perfectly conducting channels to de- 
phasing and the SNS fluctuations thus provide a sensitive probe of dephasing in a regime where 
normal transport fails to detect it. Further, the conductance fluctuations are strongly non- linear in 
bias voltage and reveal sub-gap structure. The experimental findings are qualitatively explained in 
terms of multiple Andreev reflections in chaotic quantum dots with imperfect contacts. 

PACS numbers: 73.63.Kv,74.45.+c,74.40.+k,73.23.-b 



As a consequence of the quantum mechanical interfer- 
ence of electron wavefunctions the low-temperature con- 
ductance G of mesoscopic samples fluctuates when vary- 
ing the chemical potential or an applied magnetic field. 
These conductance fluctuations were demonstrated more 
than 20 years ago as one of the first examples of meso- 
scopic quantum phenomena in sub- micron samples [1, 2]. 
Through the Landauer formula G = (2e 2 /h)J2i^i the 
conductance can be expressed in terms of a sample- 
specific set of transmission eigenvalues {%}, and the vari- 
ance of the conductance fluctuations, Var G, provides im- 
portant information about the statistical properties of 
the transmissions, such as the distribution p(T) and cor- 
relations. An important energy scale for electron interfer- 
ence in random systems is the so-called Thouless energy 
-&rh being the shift in chemical potential /i c ~ E^h suf- 
ficient to uncorrelate the transport properties. At high 
temperatures T, strong dephasing due to inelastic scat- 
tering with rate 7^ 3> i?Th subdivides the sample into 
many uncorrelated parts and the conductance fluctua- 
tions are suppressed by self-averaging. As the temper- 
ature is lowered Var G increases and it is a remarkable 
result that when T, 7^ <C -Erh (usually 7^ <C T at low 
T), Var G saturates to a value on the order of (e 2 /h) 2 , in- 
dependent of the sample size and degree of disorder. For 
this reason the phenomenon is denoted universal conduc- 
tance fluctuations (UCF) [1, 2] and in this regime trans- 
port remains practically insensitive to dephasing. 

A fundamentally different situation occurs if the leads 
to the normal (N) sample turn superconducting (S). In 
this case, a gap A opens at the Fermi level, and a sub-gap 
energy electron incident on the S interface cannot pen- 
etrate into the lead, but is instead coherently Andreev 
reflected (AR) as a hole upon injection of a Cooper pair. 
Instead of being a sum of {71}, the transport properties 
now depend on Andreev states modified by finite- voltage 
V in way highly dependent on the transmissions {7^} [3- 
71 and Landau-Zener transitions between the states lead 



to quasi-particle current [4] . These are most probable 
when levels come close for T w 1 and <fi ss 7r as schemat- 
ically illustrated in Fig. 1(a). We will show that this has 
important consequences for the statistical properties of 
the differential conductance, because its fluctuations de- 
velop extreme sensitivity to the statistics of the almost 
perfect channels, 1. 

This Letter presents the first study focused on this 
intriguing interplay of interference and Andreev pro- 
cesses and its consequences for the statistical properties 
of mesoscopic junctions. Enabled by recent progress in 
nanoscale device fabrication [8-12] we measure the low- 
temperature fluctuations of differential conductance in 
short mesoscopic SNS devices based on semiconducting 
nanowires contacted by Niobium (Nb) leads. We sys- 
tematically study the temperature and bias dependence 
of the fluctuation amplitude, the correlation potential /x c , 




FIG. 1: (a) Energy of Andreev bound states vs. phase differ- 
ence 4> f° r various T (schematic). Arrows indicate Landau- 
Zener transitions induced by the time-dependence of (f> at fi- 
nite bias, (b) dl jdV vs. V g at various temperatures (for clar- 
ity the 6.5K, 17K, 180K, and 290K traces have been off-set 
by l,...,4e 2 /h, respectively). Inset: Scanning electron mi- 
crograph of a typical device, (c) Device schematic. 
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and the average differential conductance, and find that 
the normal- lead universal limit for the fluctuations is bro- 
ken in SNS devices as was also recently pointed out by 
Doh et al. [12]. In addition we here show that unexpect- 
edly, the fluctuations maintain a strong dependence on T 
even at low temperatures (T <C i?Th) where the normal- 
state fluctuations are saturated. To explain the data we 
theoretically analyze how dephasing modifies the statis- 
tics of the almost perfect channels T w 1 and find that 
transmissions T « 1 arc suppressed. This mechanism 
explains the strong temperature dependence of the SNS 
fluctuations and shows that they provide a much more 
sensitive probe of dephasing than normal UCF. Further- 
more, varying the bias, we find that the fluctuation am- 
plitude diverges as a power-law as V — > and we ob- 
serve, for the first time, that multiple Andreev reflections 
(MAR) lead to sub-gap structure (SGS) in the fluctua- 
tion amplitude and in fj, c . The finite bias results are 
compared with computations based on MAR-theory in 
chaotic quantum dots with imperfect contacts with good 
qualitative agreement between theory and experiment. 
From this we conclude that the results are generic for 
mesoscopic SNS fluctuations. 

The nanowires are grown by molecular beam epitaxy 
and transferred to a doped Si substrate capped with 
200 nm Si02- Contacts to individual wires are defined 
by e-beam lithography, DC sputtering of 70 nm Nb fol- 
lowing a brief etch in BHF (see Refs. [10, 11] for details). 
The leads have a critical temperature T c w 1.7 K result- 
ing in a gap A = 1.76T C w 0.25 mcV at low tempera- 
ture. The inset to Fig. 1(a) shows a scanning electron 
micrograph of a typical device; the wires have diameters 
d ~ 80 — 100 nm and the distance between the contacts 
is L ~ 100 nm. The nanowires are n-type and the de- 
vice discussed here has mobility \i ~ 10 3 cm 2 /Vs, carrier 
density n~4x 10 17 cm~ 3 , mean free path l e ~ 18 nm, 
diffusion constant D ~ 60cm 2 /s, and Thouless energy 
i?Th ~ 0.4 meV estimated from the transfer character- 
istics G(V g ). Due to the design of the outer circuit, 
the measurable supercurrcnt is strongly suppressed al- 
lowing a study of the quasi-particle current alone (within 
the RCSJ/"tilted washboard" model of Josephson junc- 
tions the device constitute a strongly undcrdamped junc- 
tion). We measure the two-terminal differential con- 
ductance Q = dl/dV using standard lock-in techniques 
(V ac = 12 fiV, 77 Hz) while varying the bias V, back-gate 
potential V g (applied to the doped substrate), and tem- 
peratures from 300 K to 300 mK. In the following, data 
from one device is presented, but similar results have 
been obtained on two additional Nb-based and one Al- 
based device demonstrating the generality of the phe- 
nomena. These results and details of the device param- 
eters, the properties of the Nb contacts, and the device 
design can be found in the supplement [13]. 

Disorder in the InAs crystal together with a multi- 
faceted wire surface [14] presumably make the system 



chaotic and the barriers formed in the NS interface domi- 
nate the resistance. Therefore we compare data with pre- 
dictions from theory of MAR [4, 7] and energy indepen- 
dent scattering Random Matrix Theory (RMT) for multi- 
channel chaotic dots with imperfect contacts, see Fig. 1 
(c) [15]. This RMT is valid for both diffusive and ballis- 
tic dots if i?Th of dot + contacts is large, eV, A, T <C -Erh- 
Thus we ignore the energy dependence of % (relaxing 
this assumption makes our numerics impractical and dis- 
cussion more involved [6]). For the instructive case of 
perfect contacts we analytically find the effect of weak 
dephasing (T <C £/rh) 011 tne distribution p(T) and of 
small bias eV « A on VarCj at T = 0. For the general 
case of imperfect contacts (N = 16 channels, transparen- 
cies Tl —Tfi chosen to match the experiment) the bias- 
dependence of transport statistics is computed at T = 0. 
For details of the theory and a discussion of the role of 
contact asymmetry see Ref. [13]. 

Figure 1(b) shows examples of the measured Q(V g ) 
for V = V for various temperatures. For T < 100 K 
they exhibit a large number of reproducible, aperiodic 
fluctuations allowing a statistical analysis of the data. 
To characterize the fluctuations, we extract for each 
trace the average (Q) = (Q)v g , the variance VarCj = 
(Q 2 ) — (G) 2 , and from the correlation function F(SV g ) = 
((G(V g ) - (Q)) ■ (G(V g + 5V g ) - ((?))) the typical V g - 
scale of the fluctuations V c (proportional to \x c [16]) as 
F(V C ) = |-F(0) = ±VarCJ [17]. The normal state behav- 
ior at temperatures below T c is measured by applying a 
magnetic field B = 0.5 T to suppress the superconduc- 
tivity of the leads. 

Let us first consider the role of temperature T. Figure 
2 shows the temperature dependence of the extracted pa- 
rameters at zero bias. For T > T c = 1.7 K, (Q) is almost 
constant 3e 2 /h showing that the current is not carried by 
thermally excited carriers. At T = 1.7 K when the leads 
turn superconducting (Q) increases as a consequence of 
Andreev reflections. The increase occurs over a range 
1 K < T < T c corresponding to the T-dependence of the 
superconducting gap A(T) (included in the figure) which, 
below 1 K, is very weak and (Q) is effectively saturated. 

The fluctuation amplitude displays a different depen- 
dence on T: Upon lowering T from room temperature, 
VarC/ increases as T~ 17 (solid line). This reflects the self 
averaging discussed above and the saturation at T ~ 5 K 
agrees with _Exh ~ 5 K estimated from the transfer char- 
acteristics. The transition to superconducting leads at 
T = T c is accompanied by a sudden increase of VarCJ, 
but unexpectedly it keeps increasing all the way to the 
lowest T; the upper inset to Fig. 2 emphasizes the low- 
T behavior of Var£. Thus, the T-dependence of VarC/ 
is not governed by A(T). Interestingly for T < IK 
VarC/ seems to rejoin the T~ 17 relationship that was 
followed above 5 K. We note that the normal-state sat- 
uration value 0.09(e 2 //i) 2 (measured with B — 0.5 T) 
is of the order of the theoretical normal-state universal 
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FIG. 2: Temperature dependence of (Q) (□), VarC? (o), and 
fi c (A). For VarQ and p c solid lines are fits (aT b ). For 
VarCJ only the data points with T > 5K were included in 
the fit and the dashed extension is the extrapolation to lower 
temperatures. Solid symbols show parameters measured in 
the normal state. Upper inset shows the T-dependence of 
VarC? on log- linear scale. All data are for zero DC bias. Lower 
inset: Schematic illustration of p(T) for no (dashed), and 
weak dephasing, T, 7^ <C i?Th (shaded). 

value [18]. In the superconducting state, however, Varf? 
reaches 2.5(e 2 //i) 2 at 300 mK, ~ 30 times larger than the 
normal state value [19]. 

To describe this behavior we consider the Andreev 
states which are formed when the leads turn supercon- 
ducting. These appear at energies sensitive to the phase 
difference of the leads </> and the transparency of the chan- 
nels, e it ± = ±A(1 - 7^ sin 2 </>/2) 1 / 2 [3] as illustrated in 
Fig. 1(a). A finite bias V <C A/e leads to a quasipar- 
ticle current cx exp(— 7rA(l — %)/eV) since the result- 
ing time evolution of the phase difference, <f> = 2eVt/H, 
induces Landau-Zener transitions between low energy 
pairs e^± ~ 0[4]- Such transitions are most probable 
for 71 — ► 1 and (f> ir, and in contrast to the nor- 
mal case, transport is therefore dominated exponentially 
by the almost perfect channels. We therefore study the 
role of dephasing on the statistics of T — * 1. Using the 
dephasing-probe model [20], Ref. [21] numerically demon- 
strated that in a single-channel dot dephasing suppresses 
p(T) for T — > 1. Using this approach, we consider the 
limit 70 <C S, (S is the level spacing) and for 1 — T <C ~1<t>/& 
we have derived an exponential suppression of the trans- 
mission density, p(T) oc exp[— 0(^/26) /(l — T)] (/3 is 



the Dyson parameter). Extending to the multi-channel 
limit N 3> 1 we find that dephasing leads to the ap- 
pearance of a temperature-dependent upper bound T+ = 
1 — j^/2nE Th such that p(T) = for T > T + (see lower 
inset to Fig. 2). For normal transport this is accompa- 
nied by a practically undetectable correction — 70/71- £>rh 
to the UCF [21]. However, for SNS transport due to the 
exponential sensitivity of the current to T's near 1 the 
appearance of T + makes transport strongly temperature 
dependent and unlike normal UCF theoretically Var Q di- 
verges, VartJ — > 00, T, V — > (the F-dependence is dis- 
cussed below). In conclusion, lowering T decreases 70 (T) 
and increases T + and thus allows the exponential contri- 
butions from progressively more transparent channels to 
play a role in the transport and its fluctuations thus in- 
creasing VarC/. We are, at present, not able to predict 
the functional form of the increase and the power-law 
relationship VarC/ oc T~ 17 suggested by the experiment 
remains unexplained. Also, the combined inclusion of 
dephasing and imperfect contacts remains a challenging 
theoretical problem. 

Consider now the T-dependence of p c presented in Fig. 
2. Generally, /x c reflects the dependence of the T's on 
V g and provides information about the statistics of {%} 
complimentary to VarC?. For T > Txh the correlations 
depend on dephasing and thermal averaging and /i c de- 
creases with T and saturates for T <C Trh to a value 
~ 0.35 meV in agreement with the previous estimates of 
T-Th- As T is lowered from T c we observe a further dra- 
matic increase in the sensitivity to Vg (see Fig. 1) and p c 
decreases significantly below its normal-state saturation 
value (/x G oc T 7 ). In the S-state fi c probes correlations 
of the nearly perfect channels as discussed above. How- 
ever, the functional form of the T-dependence (and, in 
particular, the T = value of p c ) is a complicated and 
fully open theoretical problem, which needs further work. 

We now discuss the bias-dependence of the transport 
and its fluctuations. The nature of the important AR- 
processes depends strongly on V", because a sequence of 
n AR-processes that transfer a quasiparticle across the 
junction, is energetically possible only when eVn > 2A. 
At the same time, a large number of AR requires a high 
transparency. Hence, as the bias voltage is decreased an 
enhanced sensitivity to the tail p(T — > 1) is indeed ex- 
pected, and it is interesting to study the characteristics 
of the fluctuation pattern as a function of bias voltage. 
Figure 3(a) shows measurements of G(V) for two different 
V g illustrating the strong dependence on V g . Upon lower- 
ing |V| the differential conductance shows an increase at 
V ~ ±0.5 meV corresponding to enhanced quasi-particle 
transport when the peaks in the DOS of the leads line-up 
at V = ±2 A/e. The SGS at lower bias is the consequence 
of the bias thresholds for MAR as described above. 

Figure 3(b) shows a grey-scale representation of ~ 5000 
such traces covering < V g < 9.5 V. The enhanced 
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FIG. 3: (a) Q vs. V measured at 300mK for constant V g as in- 
dicated by arrows in panel (b) (Bottom curve off-set by — e 2 /h 
for clarity), (b) Grey-scale representation Q vs. V and V g 
(brighter: more conductive), (c) Experimental (Q) (□), VarC< 
(o), and \Xc (A, in units of i?Th = 5K), respectively, as a func- 
tion of \\f\. The dashed lines show the result when measured 
at 2K with the contacts in the normal state, (d) Measured 
and calculated Var<5 displayed on logarithmic scales; dashed 
line show V" 1 . Bottom curves show the residual of fitting the 
experiment to a power law oc V~°' s to enhance the SGS. (e) 
As (c) computed at T = OK. 



quasi-particle transport for \V\ < 2A/e is seen through- 
out the plot and the pattern of tilted bands of high differ- 
ential conductance observed for \V\ > 2A/e is character- 
istic of conventional gate and bias dependent fluctuations 
in disordered mesoscopic samples [22]. Extracting again 
(Q) , Var Q, and fi c for each bias value leads to the result of 
Fig. 3(c,left) which also includes normal-state data mea- 
sured at 2K. The average (Q) increases at \V\ < 2A/e 
and the SGS clearly survives the averaging. The cor- 
relation potential fi c decreases slowly as V is lowered 
from A/e, whereas Var£ displays a pronounced peak for 
V — ► 0. Both contain structure resembling MAR. To 
understand these results we theoretically consider multi- 
channel samples at T = first for eV -C A, when the 
generic distribution tail p(T -> 1) oc l/y/l-T gives a 
nonlinear dependence, I oc \J eVj A [4] . To find the fluc- 
tuations around this value in fully coherent wires or quan- 
tum dots with ideal contacts we use the known correlators 
of T,T [23] leading to Var£ oc 1/V 2 . The appearance 



of a divergence agrees with the experiment, however, as 
seen in Fig. 3(d) the experiment finds V~°' s . Therefore 
a more realistic model of the experiment should take into 
account the barriers formed in the NS interfaces. Figure 
3(c, right) shows (Q), Var£ and fi c computed for T = 0. 
fmportantly we now find VarQ close to V^ 1 qualitatively 
similar to the experimental trend. Interestingly, we see 
that the SGS for V = 2A/ne n = 1,2,... appears not 
only in the average current, but also in the fluctuations 
Var Q and fi c . The SGS peaks in /i c do not merge to form 
a divergence as for Var Q, but rather decrease slowly for 
V — > in qualitative agreement with the experiment. No 
prior theoretical results exists for fJ, c (V) and therefore the 
found agreement with the experiment is quite satisfac- 
tory. We note, that at lowest bias the computed values 
of Var Q are considerably larger than the measured val- 
ues (factor 20), and that computed MAR peaks appear 
considerably sharper and higher than observed in the ex- 
periment. We attribute this to our simplifying assump- 
tions of symmetric contacts and energy-independent elas- 
tic scattering. Most importantly, however, the strong 
suppression of Var Q induced by dephasing, is absent in 
the T = model. A quantitative agreement with the ex- 
periment is therefore not expected. Future studies could 
investigate this by repeating the measurement of Fig. 3 
at various temperatures in devices with individually tun- 
able barriers. 

In conclusion, we present the first systematic study 
of the fluctuations of SNS-transport. We find a very 
large enhancement of the fluctuation amplitude com- 
pared to normal-state UCF and an extreme temperature- 
sensitivity Var Q oc T -1 ' 7 even for temperatures where 
the normal-state fluctuations are saturated. We argue 
theoretically that this can be understood as the com- 
bined effect of the almost perfectly transmitting channels 
dominating the transport and the cut-off of transmissions 
close to one with increasing dephasing. Thus, SNS fluc- 
tuations provide a sensitive probe of quantum interfer- 
ence which might be used for measuring weak dephasing, 
unavailable from normal UCF. Moreover, we reveal that 
the statistical properties of SNS fluctuations exhibit sub- 
gap structure as a function of bias. Good qualitative 
agreement is found with numerical calculations based on 
scattering RMT and MAR theory. 

We thank J.B. Hansen for experimental support and 
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the Danish Science Research Councils (TSJ). 
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SUPPORTING INFORMATION 

Additional information relevant for the main 
manuscript is provided. Presenting experimental 
results from three additional devices, including one with 
a different superconducting metal as contact material, 
we establish our findings as general to disordered 
SNS junctions. Furthermore the details of the device 
characteristics are discussed and we present further 
details on the Random Matrix Theory used to analyze 
the finite-bias results. We discuss the construction of 
the scattering matrix of the sample and the role of 
contact asymmetry and dephasing for the statistics of 
the transmission eigenvalues. 



and details below). We note, that the gap A51 ss 115 /zeV 
of the leads for sample SI is significantly smaller than the 
Nb-based devices. Thus, at low bias the results will be 
more sensitive to thermal effects and noise as well as in- 
fluence from by the ac-bias (V ac ~ 12 /iV <~ O.lA^/e) 
required by the lock-in technique. Therefore V&rQ satu- 
rates at a larger bias <~ 0.2A for 51 than for the other 
devices. 

The figure also shows the corresponding data for sam- 
ple S4 having normal leads. As compared to S0-S3, the 
fluctuation amplitude for S4 is effectively independent of 
bias and has a magnitude similar to that of the supercon- 
ducting samples when these are biased outside the gap 
(\V\ > 2A/e). 



DATA FROM ADDITIONAL DEVICES 

In the main manuscript (MM) results were presented 
from measurements of one device (SO); an InAs nanowire 
contacted by superconducting Niobium leads. Figure 
4(a-f) shows data from one device (SI) with super- 
conducting leads based on a Ti/Al/Ti trilayer[ll] with 
A w 115 /xeV and two additional Niobium-based devices 
(S2,S3). The qualitative behavior analyzed in relation to 
Fig. M3 (in the following references to figures in the main 
manuscript is written as Fig. Mx) is again observed for 
all three devices with the key features being: 1) The en- 
hancement of (Q) at the quasiparticle onset \V\ — 2A/e, 
and sub-gap structure at lower bias. 2) The strongly 
peaked fluctuation amplitude as |V| — > OV (see below) 
and structure in Var Q at lower bias, and 3) The slow de- 
crease of the correlation potential V c for \V\ — > OV and 
sub-gap structure also in V c . Furthermore, for compari- 
son, Fig. 4(g-h) shows the corresponding measurements 
for a similar nanowire device contacted by normal-metal 
Titanium/Gold leads[10]. The low-bias behavior of (Q), 
Var Q, and V c remains featureless, thus confirming the 
significant role of the superconductors. The increased 
noise of the results in Fig. 4 in comparison with the re- 
sults of the main manuscript is due to poorer statistics 
for S1-S3. Nevertheless the results supports the generic 
nature of the discussed phenomena. 

In Fig. 5 the bias dependence of the fluctuation am- 
plitude for samples S0-S4 is displayed on logarithmic 
scales where each data point represents an average of 
positive and negative bias values. For the samples with 
S-leads (S0-S3) the Var£ oc V v dependence discussed 
in the main manuscript is again observed with best-fit 
powers —1.13 < p < —0.79 (inset) in good agreement 
with our numerical results for the fluctuations of multi- 
channel quantum dots with imperfects contacts (details 
below). We ascribe the spread in p to sample-specific 
transparency and asymmetry of the contacts which have 
great influence on the fluctuation amplitude and its bias- 
dependence (see the discussion in the main manuscript 



DEVICE PARAMETERS 

The characteristic parameters of the devices (carrier 
density n, Fermi wave vector kp, Fermi velocity vf, mean 
free path l e , number of channels N, diffusion constant D, 
and Thouless energy i?Th) were estimated from the mea- 
surements of the transfer characteristics and bias spec- 
troscopy in the Coulomb blockade regime. In the fol- 
lowing we describe the analysis for sample 5*0 but the 
parameters of 51 — 54 are obtained similarly and the 
values arc collected in Table I. Figure 6(a) shows the 
linear conductance G vs. V g at a temperature of 17K 
where conductance fluctuations are not yet dominant 
(The analysis is insensitive to the temperature as the 
G vs. V g traces for T < 150 K are very similar. Below 
~ 10 K, however, UCF makes an accurate determina- 
tion of the transconductance problematic). As seen in 
Fig. 6(a) the wire is depleted from carriers at low gate- 
voltages. The threshold appears at Vq,t ~ — 8 V from 
which on G increases linearly with V g until V g ~ —5 V. 
Within the charge control modcl[24] , which is widely used 
for analysis of nanowire FET's[25], the transconductance 
g m = dG/dVg w 0.6e 2 /hV is given by g m = [iC g /L 2 7 
where /i is the mobility, C g the capacitance to the back- 
gate, and L <~ 100 nm the length of the device. The ca- 
pacitance is found from the V^-separation SV g ~ e/C 9 of 
Coulomb blockade conductance peaks which appear close 
to pinch-off[10] (this slightly underestimates C g due to a 
finite level-spacing, however, this correction is not sig- 
nificant for the analysis). Figure 6(b) shows a measure- 
ment of dl I dV vs. bias and gate for V g ~ — 7 V exhibit- 
ing the characteristic Coulomb-blockade diamonds with 
SV g ~ 65meV yielding C g s=s 2.5 aF in good agreement 
with similar studies in other nanowire devices [10]. This 
C 9 -value is somewhat lower than the result of an electro- 
static cylinder-over-plane model where C g — ^(ihjr) ~ 
9aF (eo,e r ,L,r,h are the free-space permittivity, relative 
dielectric constant of SiC>2, length of wire-segment be- 
tween the leads, the wire radius and center-to-plane dis- 
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FIG. 4: Measurements from additional devices: (a),(b) Results from a nanowire device (SI) with a Ti/Al/Ti tri-layer contact 
turning superconducting below 750mK (A « 115 ^teV). (c),(d) and (e),(f) Results from nanowire devices (S2 and S3) with Nb 
contacts similar to those of the main manuscript (SO). (a),(c),(e) Conductance as a function of bias and gate voltages showing 
the increased conductance below the superconducting gap and the characteristic pattern for gate and bias dependent CF. 
(b),(d),(f) Average conductance, variance, and critical voltage as a function of bias voltage, displaying similar phenomenology 
as displayed in Fig. M2 for device SO. (g),(h) Results from a nanowire device (S4) with normal Ti/Au contacts to allow a 
comparison with the corresponding measurement in the absence of superconductivity. All results are measured at T = 300 mK. 



Dev. p 




t < < < — ■— 

0.1 

V[A/e] 



FIG. 5: Variance vs. source-drain bias V in units of A/e 
for samples S0-S4 on logarithmic scales. Positive and nega- 
tive bias values have been averaged. A power-law dependence 
VarC/ oc V p is observed for devices S0-S3 which have super- 
conducting leads and the best-fit exponents are collected in 
the table. In contrast, the results of S4, which has normal- 
state (Ti/Au) leads display a strong saturation (for S4, A was 
set to 0.25 meV to allow easy comparison with the Nb-based 
devices S0,S2, and S3). All results are measured at 300 mK. 



tance, respectively). However, this is expected since the 
electric field from the back-gate in our device is effec- 
tively screened by the large metal leads (see SEM-image 
on Fig. Ml). As C g depends mainly on the geometry of 
the sample we do not expect it to depend significantly on 
V g . Using C g — 2.5 aF we get a mobility \i ss 960 cm 2 /Vs. 
The carrier density in the wire at a gate potential V g can 
be estimated from the charge induced by the gate- voltage 
with respect to the pinch-off, n = C g (V g ~ Vg ,t) / nr 2 Le 
(here r is the wire radius and L the length of the wire 
segment between the leads). This gives n « 2.5 — 5.7 x 
10 17 cm~ 3 for Vg = - 10 V. We note that these esti- 
mates of mobility and density are similar to other studies 
of InAs nanowire devices[8, 26]. 

The Fermi wave vector Uf is calculated using the 3D ex- 
pression for the Fermi energy fcp- = (6n7r 2 ) 1 / 3 and using 
the bulk value for the effective electron mass in InAs 
to* = 0.026m e (m e being the electron mass) we find 
also the Fermi velocity vf — http/m*. Finally, since 
fi = ele/vFTn* we find the mean free path l e , the dif- 
fusion constant D = -gVple, an d the Thouless energy 
Et\i = hD/L 2 . In order to transform changes in V g into 
the corresponding change in chemical potential the pa- 
rameter a = Cg/Ctotai is needed. In analogy with the 
standard procedure for finding a in the Coulomb block- 
ade regime it is determined as the typical slope of the 
high-conductance ridges observed in Fig. M3(b) giving 
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(a) 




0-6 e ? /hV 




T=17K 





(b) 
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FIG. 6: (a) The linear conductance as a function of V g mea- 
sured at 17 K for device SO. For V g above the threshold 
at Vg,t ~ — 8V the conductance increases linearly with a 
transconductance g m « 0.6e 2 //i1/ which allows a determina- 
tion of the mobility, (b) as (a) but measured at T = 300 mK 
and for V g close to pinch off where Coulomb Blockade domi- 
nate the transport. 
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FIG. 7: (a), (b) Grey-scale representation of dl/dV vs 
and V g for applied magnetic field B — T and B = T, 
respectively (brighter, more conductive), (c) dl /dV vs. V and 
temperature at fixed gate potential. The BCS-like evolution 
of 2 A is indicated by the dashed line, (d) Traces of dl/dV 
vs. V for various combinations of temperature and magnetic 
field as indicated on the figure. 



a = 2.5meV/V for SO. The results are summarized in 
Table I which includes also the results of a similar anal- 
ysis for devices S1-S4. 



to 2 x 1.752c = 0.5 meV at low temperature. This 
agrees with the observed increase of the conductance for 
|V| ~ 0.5 mV as the peaks in the density of states of the 
leads line up at 2A/e. 

The critical temperature measured in the device is 
considerably lower than that of bulk Nb which has a 
critical temperature of 9.2 K. Such differences between 
bulk properties of the lead material and the actual 
properties of the nanoscale devices is often observed. 
For example, in Refs. [11, 27] aluminum was used 
for contacting carbon nanotubes and InAs nanowires, 
respectively, with an observed transition temperature of 
750 mK considerably lower than T c of bulk Al (1.2 K). 
In the case of the present device the dramatic decrease 
of T c may be due an impure Nb film in the interface 
between contact and nanowire due to a reaction of the 
sputtered Nb with outgassing from the electron-beam 
resist (PMMA) on the substrate or traces of oxygen 
during sputtering[28, 29]. 

The critical magnetic field of Niobium can exceed several 
Tesla depending on the quality and geometry of the film, 
making Niobium a good candidate for nano-structure 
based SQUID's[30] where robustness to an external 
magnetic field is desired. In our case, however, the 
critical field turns out to be relatively small which is 
beneficial as we can then measure the normal-state 
behavior for temperatures below T c . 
Figure 7(a) and (b) shows dl/dV vs. V and V g with 
B = 0T and B = 0.5 T, respectively. For B — the 
conductance exhibits pronounced sub-gap peaks due to 
Andrecv reflections (similar to Fig. M3(b)). The sub-gap 
structure disappears upon the application of ~ 350 mT, 
however, at B = 0.5 T, a small conductance decrease 
remains for |V| < 0.1 mV as also seen in Fig. 7(d). This 
feature repeats for all gate-voltages and indicates that 
some reminiscence of superconductivity may still exist. 
From the measurements for T > T c in Fig. M3 and 
of sample S4 in Fig. 4 it is known that in the normal 
state the statistical properties (Var (?,((/) and fi c ) are 
bias-independent for small V. Therefore, to ensure that 
the normal-state data reported for T < T c in Fig. M2 
(solid points) are free from superconducting correlations 
they are measured with B = 0.5 T and a small bias 
V = 0.2 mV. 



PROPERTIES OF THE NB FILM 



The results reported in the main manuscript are from 
a device with an InAs nanowire contacted by Niobium 
leads deposited by DC sputtering. In Fig. 7(c) shows 
traces of dl/dV vs. V at constant gate V g = 2.5 V as 
the temperature is lowered from 1.8 K to 0.3 K. Below 
T c ~ 1.7 K the low-bias conductance increases and 
at lower temperatures discernable sub-gap structure 
develops. The dashed curve shows the temperature 
dependence of (twice) the gap 2Abcs(T) saturating 



CONSIDERING THE SUPPRESSED 
SUPERCURRENT 

The theoretical supercurrent through the wire is 
given by [31] I c = nABCs/^eR n where R n is the 
normal state resistance of the device. In our case 
R n ~ (3e 2 //i) _1 which yields a theoretical supercurrent 
of I c ~ 30 nA. The measurable supercurrent, however, 
depends on the external circuit and in the extended 
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Sample id. 




SO 


SI 


S2 


S3 


S4 


Lead material 




Nb 


Al 


Nb 


Nb 


An 


Device length [nm] 


L 


100 


300 


100 


100 


300 


Device diameter [nm] 


d 


80 


70 


80 


80 


70 


Gate capacitance [aF] 


c g 


2.5 


1.2 


2.5 


1.6 


1.8 


Carrier density [10 17 cm -3 ] 


n 


4.1 


1.9 


4.7 


2.0 


1.2 


Mobility [10 3 cm 2 /Vs] 


/' 


1.0 


4.1 


0.9 


1.0 


4.5 


Mean free path [nm] 


le 


18 


60 


18 


15 


56 


Fermi wavevector [10 6 cm -1 ] 


kF 


2.9 


2.2 


3.0 


2.3 


1.8 


Fermi velocity [10 8 cm/s] 


Vf 


1.3 


1.0 


1.3 


1.0 


0.8 


Fermi wavelength [nm] 


Xf 


22 


28 


21 


27 


33 


Diffusion constant [cm 2 /s] 


D 


60 


180 


73 


45 


120 


Thouless energy [meV] 


-ETh 


0.4 


0.1 


0.5 


0.3 


0.1 



TABLE I: Results of the analysis based on the transfer characteristics of the nanowire device. The analysis is based on the 
charge control model and the quoted values are the mean values over the relevant gate voltage interval. 



RCS J/" tilted- washboard" model[32] it depends on the 
quality factor Q, where Q^ 1 = uu p (RC + I( } R ) 

and lu p = ^2eI c /h{C{l + R/R n ) + Cj). Here Cj is 
the junction capacitance and Ij = I c sin(0) is the 
phase-dependent super-current through the junction. 
C is the relatively large area bonding-pad capacitance 
and R is the resistance of the on-chip wiring connect- 
ing the bonding pads to the device. The circuit is 
schematically shown in Fig. 8, where Rw ^ 50 £1 is the 
resistance of the wiring of the cryostat. In Ref. [27] 
efforts were made to optimize these parameters for a 
large measurable supercurrent in carbon nanotube based 
junctions (essentially by making R large). In the present 
case of the Nb-NW-Nb junctions we have Cj ~ 2 £F, 
R n ~ 10m, C ~ 1 - 2 pf , and R < 10ft, and with these 
parameters we get Q ~ 10 showing that the measurable 
supercurrent is strongly suppressed. This is consistent 
with the experimental results: At all gate voltages 
we observe the reminiscence of the highly suppressed 
supercurrent as a narrow weak peak in dl/dV at zero 
bias (see Fig. M3(a)). The peak is, however, negligible 
compared to the quasiparticlc conductance which allows 
us to analyze the results in terms of the quasiparticlc 
transport and Andreev reflections alone. 



POSSIBLE EFFECTS OF THERMAL NOISE AND 

Vac 



cryostat - 

-- on-chip 




FIG. 8: Schematic circuit-diagram of the device and measure- 
ment setup. Cj, R n , Ij: Junction capacitance and normal- 
state resistance and current, respectively. C is the bonding- 
pad capacitance and R is the resistance of the on-chip wiring 
connecting the bonding pads to the device. 



Vat(0.3K) < 12 fieV, so the Nyquist contribution must 
satisfy V N (T) < y/T[K] x 22 ^eV. 

The approximate form of Var Q as a function of voltage, 
VarG/{e 2 /h) 2 w 0.23(A/eU)°- 8 for V > V ac ,V N (0.3), al- 
lows us to evaluate the upper boundaries for VavQ al- 
lowed in our experiment at V = 0. Using V = V ac 
and A = 250 [ieV gives 2.6, while V = V N (T) gives 
> 1.6T~ a4 ,T > 0.3 K. Only at lowest temperature do 
our data reach 2.6, and at higher T values significantly 
lower than the Nyquist upper bound are observed. Thus, 
we conclude that the power-law extracted from our data 
cannot be governed by a thermal noise. Further, the 
values of Var (J are affected by V ac only at the lowest 
temperature. 



In principle the temperature dependence presented in 
Fig. M2 could be affected by thermal noise and the finite 
excitation voltage of the lock-in amplifier V ac = 12/j.eV. 
Due to down-mixing by the nonlinear SNS device, ther- 
mal noise may result in a dc voltage Vjv(T) oc VT. 
Since a bias-dependence of Vart/ is observed in Fig. 
M3 at T = 0.3 K down to V ~ V ac , we conclude that 



CONSTRUCTION OF 5-MATRIX 

Since we expect that the fluctuational phenomena we 
consider are generic, particular geometry of the sam- 
ple is not expected to make qualitative difference. On 
one hand, the barriers on the NS boundaries are the re- 
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gions where the main voltage drop occurs. The meso- 
scopic sample itself can then be considered point-like. On 
the other, the samples are disordered due to impurities. 
These two features allow us to consider a chaotic quan- 
tum dot as a realistic geometry for the experiment. Im- 
portantly, the quantum dot model allows us to study how 
the fluctuations depend on the voltage V g on a nearby 
electrostatic gate and on the imperfect contacts. Such 
information, unavailable from other models, is important 
for the comparison to the experiment: 

1) Effect of the gate voltage V g . Conventionally, a 
closed chaotic dot is characterized by its large M x M 
Hamiltonian matrix TL from a Gaussian Ensemble of ran- 
dom matrices with relevant symmetry. This symmetry is 
characterized by Dyson parameter for pure ensembles: 
Orthogonal, /3 = 1, if time- and spin-reversal symme- 
try are present (GOE), Unitary, = 2, if time-reversal 
symmetry is broken (GUE), or Symplectic, /3 = 4, if 
spin-symmetry is broken (GSE). The ensemble average 
(denoted by (...)) of any Hamiltonian element vanishes, 
(Ha-y) = 0. The pair correlator of the elements is defined 
by /3, the matrix size M, and mean level spacing <5, [33] 



contacts, see Fig. 9: 



7 'I-q / 7 / / 



9 {^a-f'S^a' + 0~aa'Sy~fiO~l3l) , /? — 1, 2. 

7T 



An open dot with N ballistic channels is fully charac- 
terized by its N x N scattering matrix U. Usually hi is 
assumed to be uniformly distributed over the ensemble 
of unitary matrices of relevant symmetry (Dyson Circu- 
lar Ensembles with (3 = 1,2,4). However, a gate voltage 
V g coupled to the dot via capacitor C, see Fig. 9, shifts 
the chemical potential (the bottom of the band) and thus 
affects the matrix U. Variations in V g lead to a shift of 
the dot's Hamiltonian and therefore 



U 



t N - 2 7 rilU t 



W. (1) 



The unit matrix of N x N size is denoted by Iat and 
the coupling M x N matrix W consists of the matrix 
(y/M5 /tt)1n and zeros in the lower M — N rows [33]. 
For N ^> 1 universal results (independent of M) are 
reached only when M — + oo. In the limit N -C M the 
Hamiltonian Gaussian Ensemble was shown to give the 
same uniform distribution for IA as the Circular Ensem- 
bles [34]. Equation (1) is essential for finding the role of 
V g in numerical mesoscopic averaging. 

2) Effect of imperfect contacts. For ballistic contacts, 
the scattering matrix hi is taken from the Circular En- 
semble directly, or by using Eq. (1). For a dot with im- 
perfect channels with transparency {^},0 < G% < l,i = 
1,..,N the scattering matrix S is distributed according 
to the Poissonian ensemble [15]. A representative of this 
ensemble can be obtained from the matrix U of an open 
dot after including possible multiple reflections from the 



s = y/i N -g -Vgu 



1 



In - V^N - QU 



Vg. (2) 



Indeed, expanding the last term to n-th order in U ac- 
counts for ri — 1 internal reflections before the electron 
exits through one of the contacts. 

We use quantum dots to model the SNS samples and 
find the fluctuations of conductance G = dl/dV. Each 
sample is specified by its (S-matrix, and we find its set of 
transmission eigenvalues {%}■ Using the scattering the- 
ory which includes multiple Andreev reflections (MAR) 
theory in SNS structures, developed by Averin and Bar- 
das [35] , we compute the current / as a sum of currents in 
each channel for the slightly shifted voltages to find the 
sample-specific G. Repeating this calculation for many 
different 5-matrices allows to find statistical properties 
of G. 

In general, the mesoscopic averaging can be performed 
after measurements in many samples, or by using a sin- 
gle sample and varying the gate voltage V g . Justification 
for this widely used procedure comes from the hypothesis 
that energy averaging equals ensemble averaging. There- 
fore, the mesoscopic averaging denoted by (...) should be 
understood as (—)v' g - 

Experimentally, correlations for traces taken at the 
same bias voltage V and gate voltages shifted by 5V g 
are quantified by the correlator F [17], 



F(V,SVg) = (G(V,Vg)G(V,Vg + SVg)) Vg 

-(G(V,V g )) Vg (G(V,V g + 6V g )) Vg 



(3) 



The variance of the conductance is then given by Var g — 
F(V, 0). For large SV g the conductances become com- 
pletely uncorrelated, F(V,6V g ) — + 0, and the correlation 




FIG. 9: Imperfect contacts for a dot with ./Vi,2 channels are 
modeled by transmissions Q\ t 2 for the contacts. The S ma- 
trix in Eq. (2) accounts for multiple reflections due to <5 li2 7^ 1. 
The gate voltage V g varies the potential in the dot via a ca- 
pacitor C g . The coupling of the dot to a probe with N$ chan- 
nels and transmission Q$ corresponds to the dephasing rate 
70 = N^Q^S. Uniform dephasing for a given 7^ is reached at 
iV^, — » 00, — » 0. 
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potential p c is defined as the shift in chemical potential 
jj, c which diminishes the correlator twice, 



F{V,6V g °) = ±F{V,0), Mc 



aSV n c 



(4) 



In normal transport through a dot with ballistic con- 
tacts ji c is naturally measured in units of its escape 
rate, or the Thouless energy of an open dot, _E Th = 
NS/2tt — h/rd- For example, for normal linear conduc- 
tance p c = i?Th , = 1,2 and it reaches its maximum 
[(2V2 + 3i/2)( 2 V2_i)]V2£ Th ~ \,UE T h in the crossover 
between = 1, 2. If the contacts are imperfect with equal 
transmission Q, the Thouless energy (and thus p c ) dimin- 
ishes according to the time an electron typically spends 
inside the dot, E Th = NGS/2tt. 

In SNS transport, to find p c (V) we have to numerically 
solve Eq. (4) (iteratively by the Newtonian method) . We 
substitute the gate voltage averaging in Eq. (3) by averag- 
ing over several hundred (100-400) samples at each iter- 
ation, their S matrices being found combining the proce- 
dure in l)-2). While the correlator F is mathematically 
well-defined in Eq. (3), in rc£ dity both sides of Eq. (4) 
fluctuate, and our procedure may converge slowly. Con- 
vergence after few iterations results in noise in V&rQ(V), 
see Figs. 3(d,e) in the main text. Indeed, averaging the 
r.h.s. of Eq. (4) in just few hundred samples may differ 
from the true value of V&rQ. Similar fluctuations in the 
l.h.s. result in noise for the p c (V) plots. 

The Random Matrix Theory (RMT) described here 
uses energy-independent S matrices for electrons with 
(generally) different kinetic energies. This assumption 
is valid if the electrons are close to the Fermi level. A 
typical energy of an electron is limited by the largest 
energy among the temperature scale ksT, the supercon- 
ducting gap in the contacts, A(T), and the bias eV. The 
scattering matrix S can be taken energy-independent if 
T, A, eV <C ETh — NQ5/2ir due to large level spacing 6 
and good conduction of the sample, NQ 3> 1. 

In addition, for validity of our analytical results wc 
must ensure that for small bias voltages, eV <C A, 
the band <~ (eV/A) 1 / 2 still contains many transmission 
eigenvalues T. This is fulfilled, if 1/N < (eV/A) 1 ' 2 < 1. 
Numerics for N 3> 1 are performed for = 2, but the 
results are easily generalized on — 1,4, since p(T) is 
insensitive to for N » 1 and the correlator JC is simply 
rescaled with 1/0 [23]. 



EFFECT OF IMPERFECT CONTACTS 

Analytical results for Var£ in SNS transport can 
be obtained after combining transmission correlators 
JC(T,T') [23] with the MAR theory of Avcrin and Bar- 
das [35], who showed that for eV « A only channels 
with high transmissions, T w 1, are important. Landau- 
Zcncr transitions between Andreev bound states lead to 



nonlinear current / oc (eV/A) 1/>2 in diffusive wires [36]. 
We point out that p(T) oc 1/ yJl—T is generic for meso- 
scopic samples with perfect connection to the leads, such 
as diffusive wires or quantum dots with ballistic contacts, 
and obtain 



(I(V)I(V')) ex 



Wv 7 

0(V + V>) 



(5) 



Var Q is obtained from Eq. (5) after differentiating both 
sides of equation dydyi and setting V' = V. The fi- 
nal result is V&rQ oc l/V 2 ,eV <C A. Numerical results 
for multi-channel dots with ballistic contacts, N = 16, 
indeed show this behavior. However, this instructive ex- 
ample does not take the contacts into account, which are 
important for the current experiment. 

The opposite limit is a random mode mixer, or the 
Fabry-Pcrrot interferometer, where electrons gain ran- 
dom phases traversing between the contacts. The in- 
ternal reflections are absent and the main resistance 
comes from the contacts ('opaque mirrors') [37] [38]. This 
model is thus relevant for almost perfect conductors and 
its results are also universal. The distribution p(T) of 
such a random Fabry-Perrot interferometer depends on 
the transparency of the contacts, and is bound by 

r_ <r <r+ [37], 



P{T) 



G1G2 

(1 T — £i)(l — G2)) 2 



(6) 
(7) 



The transmission can be parameterized by a uni- 
formly distributed phase ip G [0, 2tt] as T = 
T + T_/(T + cos 2 ip/2 + T_ sin 2 ip/2). Even though the cor- 
relator K,(T,T') oc 1/(0 sin 2 {ip — (f')/2) is formally inde- 
pendent of the cut-offs 7± , they do affect fluctuations for 
eV <C A: 



(I(V)I(V')) 



7T0(V + V) 

'nA(T+-l)(V + V') 



x cxp 



eVV 



■ (8) 



Asymmetry in transmissions of the contacts Q\ ^ Q 2 
modulates 7± in Eq. (7) and suppresses currents and their 
fluctuations. Importantly, the appearance of a fixed cut- 
off for perfect channels affects the current fluctuations 
exponentially. One reason for the appearance of 7+ < 1 , 
the contact asymmetry, is obvious from the last exam- 
ple, and below we consider it more quantitatively for 
chaotic quantum dots. The universality of the results 
(5), (8) does not hold in a general situation: the shape 
of p(T) and the cut-off values 7±, found by the meth- 
ods of Refs. [39], [18] depend on conductance and the 
contact asymmetry, N\ =/= N 2 or Q\ ^ Q 2 . Symmetric 
contacts with N x = N 2 still yield p(T -> 1) oc 1/y/l-T, 
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even for Q\ = Q 2 < 1, 



but for asymmetric contacts we 
For transmission 



can have p(T -> T+ - 0) oc - T. 
distribution p(T) we can take into account the contact 
asymmetry following the method of Rcf. [18]. We find 
that the perfect channels do not vanish, p(l) ^ 0, if the 
contacts are not very asymmetric, and p(T) near T = 1 
is suppressed compared to the universal distribution as 



P{T)\ 



^]{Ng 1 -gN 2 ){Ng 2 -gN 1 ) 



T- 



irgVT^r 



(9) 



where Q = ^1+^2-^1^2- Indeed, for N ia = N/2 perfect 
channels exist when the contact transparencies are close, 
\Q\ — Q2\ < g^2- In the limit £i, 2 <C 1 this condition can 
be easily violated by a relatively small difference in g. We 
can not analytically predict the behavior of the current 
fluctuations, since /C(T, T') is unknown and most proba- 
bly non-universal. However, we assume that our contacts 
are not very asymmetric and p(l) 7^ 0. Even if asymme- 
try in contacts does affect Var g in our experiment, it can 
not account for a strong temperature dependence of our 
data. To explain the strong T-dependence of our results, 
we consider the effect of dephasing on the transmission 
statistics.. 



EFFECT OF DEPHASING 

As discusses above, in the low-temperature limit, 
70, T <C ETh, and small bias voltages, eV <C A, the 
statistics of T close to perfect transmission become ex- 
tremely important for SNS transport. As seen in Eq. (8), 
the appearance of a cut-off T+ < 1 strongly suppresses 
VarC/. Our samples are expected to combine effects of 
imperfect contacts and dephasing. For simplicity we now 
consider the role of dephasing alone, assuming ballistic 
contacts to the reservoirs. Temperature is assumed to 
be sufficiently low to make the dephasing rate small, 
70 = h/T<f, <C ETh- The dephasing probe model, pro- 
posed by Biittikcr [20] , has been extensively used in the 
literature. In this model the quantum coherence in the 
dot is destroyed by attaching a probe, where the volt- 
age can be either externally controlled or remain float- 
ing. This probe exchanges electrons with the dot, and 
one distinguishes non- uniform (or localized) vs. uniform 
dephasing by such a probe depending on its coupling to 
the dot, see Fig. 9. 

Uniform dephasing denotes a probe having large num- 
ber of poorly coupled channels, Q^, — > 0, — > 00. How- 
ever, the dephasing rate 70 = N^g^S remains fixed. 
Known results for uniform dephasing correspond to the 
results of the imaginary potential model in the Hamilto- 
nian approach [21]. Non-uniform dephasing, on the other 
hand, usually takes perfect coupling to the dot, = 1, 
and finite N^. The advantage of this approach is the 
technical simplicity due to lack of back-reflection from 



the probe, but on the other hand 70 should be strictly 
quantized. For dots which are not artificially dephased 
by a local probe we assume the uniform dephasing model 
to be more realistic. 

For small dephasing in a dot with arbitrary TV only the 
uniform dephasing model can be used. To illustrate the 
role of 7^ in the tails of p(T) we first take a single-channel 
quantum dot, N — 2, and later consider the limit JV>1. 
Coherent quantum dots with N = 2 have only one trans- 
mission eigenvalue T and p(T) = (f3/2)T^ 2 ~\ (3 = 1,2 
[23]. For weakly dephased dots we use the intermediate 
results of Ref. [21] and express the dimensionless con- 
ductance g = hG/2e 2 close to g = 1 as a sum 



2 x ^ x j 



M 1 - 

1,3 = 1 ■> 



u u u i2 u 1:j u j2 



(10) 



where u, u' are 2x2 unitary random matrices and u' = 
u T , (3 = 1. If the particles were absorbed by the probe, 
only the first term in the sum (10) would have been 
present. The second term is due to reinjection of par- 
ticles by the probe and results from the requirement of 
particle conservation. The parameters N^Q^/xi^ char- 
acterize the coupling strength of the probe to the dot, 
x — > 00 corresponds to weak coupling. The distribution 
p(T) for N^g^ <C 1 is found after integration over the 
uniform distribution of u, u' in the unitary group and the 




2JV(1-T) 



FIG. 10: The tail of the distribution p(T) in a weakly de- 
phased quantum dot, N^Q^ <C 1. Distributions in a single- 
channel dot, N = 2, plotted for j3 = 1 (long-dashed) and 
(3 = 2 (solid) are exponentially suppressed at 1 — T -C NtpGj, 
and return to p ~ /3/2 at 1 — T ~ N^Q^. They should be con- 
trasted with qualitatively plotted p(T) oc (T + -T) 1/2 , T k T+ 
(short-dashed) for a multi-channel dot, N 3> 1. 
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distribution P(xi,x 2 ) (see Ref. [21] for its general form), 
p(T) = / dx 1 dx 2 / dudu'S(T — g)P(x\,x 2 ), 

p( XUX2 ) « fl*i*'i*i-*»i)y **i + « a )/2. (ii) 



48 

For 1 — T <C iV^C/^ we find an exponential suppression, 
p{T) oc exp[-/37V £ /2(l - T)]. This result p(T) cx 
cxp[— {3W(T)] can be understood as the density of clas- 
sical particles at a point T e [0, 1] and temperature 1//3 
in an external potential 



W(T) = 



2(1 -T)' 



(12) 



Even though this potential is weak, it repels the trans- 
mission T from the perfect T = 1. Only at iV^C?^ <C 
1 - T « 1 does the density return to ~ 1, see Fig. 10. 
The distributions reach half of p(l) — 0/2 at 1 — T w 
0.17A^, /3 = 1 and 0.20iV^,/3 = 2. This result can 
be interpreted as an effective cut-off T + in transmissions 
due to finite dephasing, 1 — T + <~ N^Q^ — j^/S. 

Can we expect that in multi-channel limit, iV 3> 1, 
the density p(T — s- 1) is enhanced in dephased quantum 
dots? This question is natural from comparison between 
p(l) = 0/2, N = 2 and p(T) oc 1/ A /T _ T -» oo,iV > 1 
for T — > 1 in coherent dots. It turns out to be more con- 
venient to consider not T e [0,1], but A > 0, such that 
T = 1/(A+1). Following our result for N — 2, we assume 
that some potential W(X) induced by dephasing acts on 
each eigenvalue A. The exact form of this potential and 
its dependence on dephasing strength are unknown. De- 
phasing potential W(X) should be distinguished from a 
many-body potential Wo (A) and repulsion w(A, A') from 
another eigenvalue A'. The former, W — (N/2) ln(A+l), 
weakly repels A from oo and appears due to smaller phase 
space available to large A in the limit N ^> 1 (we ne- 
glect with corrections 0(1) to Wo(A), leading e.g. to 
weak localization correction). The latter, the repulsion 
u(A, A') = — In |A — A'| between the eigenvalues A, A', ori- 
gins in the chaotic dynamics in the sample, and we fur- 
ther assume for simplicity that it maintains its universal 
form [23] (see discussion below). 

If W(X) + Wo(A) and/or asymmetry in the contacts, 
Ni 7^ N 2 , lead to cut-off values A = a, b such that p(X) = 
0, A = a, 6, we find [40] 



P(A) 



1 



Nny/(a + l)(b+l) 



ny/(X-a)(b-X) 



A + l 



-/ 

J a 



b dX ,dW(X')^ 



a){b-X') 



dX' 



A - A' 



+ 7TZ 



(13) 



where z = N m i D —N/2 < is defined by the contact asym- 
metry. The distribution p(X) depends on the exact shape 
of W(X). Integration of (13) over A satisfies the nor- 
malization condition, J dXp(X) = iV m i n , since the W(X') 



contribution is canceled after integration over A. When 
W = 0, the boundary is given by V&o + 1 = —N/2z, 
and the minimal possible value of A, an = 0. For sym- 
metric dots, z — 0, one reproduces bimodal distribu- 
tion p(T) = N/(2y/T{l - T)). However, if the poten- 
tial W(X) is unknown, the resulting distribution p(X) 
must reproduce results obtained in the uniform dephas- 
ing model, VarCJ in particular [21]. The transmission re- 
pulsion of coherent dots allows us to add another equa- 
tion relating VarCJ and a, b: 



Varq 
(2e 2 /h) 2 



T) 2 = 



^N 2 



1 



a+l 

2 



6+1 



N(N + N^g^) 



(14) 



Together with boundary conditions on p(X) the total sys- 
tem of equations reads 



I" 

I" 

J a 



dW(\) 



b dMX+1) dW{\) 
V(A-a)(b-A) dx 

_L 1_ 

1+a 1+6 



2^/(l+a)(l+6)' 
TTZ, 
N 2 -4z 2 



(15) 



If the potential only repels A from 0, the second and 
third equations in (15) for symmetric dots, z — 0, give 
a = N^g^/Njb = oo. If we additionally assume that 
W(X) oc N^Qcf,, the first equation in (15) results in 
W(X) — A^CJ^, arctan(l/v / A)/"\/A. However, a differ- 
ent assumption about functional dependence of W(X) on 
N^Qff, leads to different dependence on A. Therefore, ad- 
ditional information about IF(A) is needed. 

On the other hand, we could assume that the poten- 
tial W(X) maintains the functional form given by (12), 
W(X) = (A^/2)(l + 1/A). This potential only repels A 
from the perfect transmission, A = 0, and for a symmet- 
ric dot, z = 0, the second equation in (15) gives b — > oo. 
First equation in (15) allows us to express a as a function 
of dephasing strength and find the total distribution: 



P(A) 



a + 1 



AVX~^[(1 + 2a)A + 2a] 
27r v / aTT(A + 1)A 2 



(16) 
(17) 



Instead of generic p(T) oc 1/yT — T,T — ► 1 the distribu- 
tion behaves as p(T) oc (7+ -T) 1/2 ,T «T + , see Fig. 10. 
However, the cut-off value a in Eq. (17) does not satisfy 
the third equation in (15). In fact, Eq. (12) is only an 
asymptote A — > and the part which repels A from oo 
is absent. This part might become important for strong 
dephasing, see numerical results of Ref. [21]. The poten- 
tial (12) leads to formation of the gap in transmissions, 
p(T) = 0, T + < T < 1, but it probably overestimates the 
size of this gap. 
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The system of equations (15) defines a, b through the 
unknown potential W(X). The function W{\) cannot be 
universal, since in the strong dephasing limit N^Q^ » N 
the last equation gives bsjI, but then the first two equa- 
tions can not be simultaneously satisfied. This shows 
that, in general, the phenomenological treatment of de- 
phasing in form of a wisely chosen potential TK(A) is not 
complete. The repulsion between transmissions should 
itself change from its universal form — In | A — A' | as a re- 
sult of dephasing and a more accurate treatment should 
be taken. 

However, from the analytical results for N = 2 we ex- 
pect that the potential W(X) is the main effect of weak 
dephasing, N^Q^ <C N and modification of transmission 
repulsion would affect the size of the gap only pertur- 
batively with small parameter N^Q^/N -C 1. Without 
finding the exact form of W(X), we find from the third 
equation in (15) perturbations to the coherent values of 
the boundaries ao,6o- If the contact asymmetry affects 
only 71 , we find 

N 4z 2 
T+ ~ N + N+G,/,' T - ~ N(N + N^Y (18) 

Comparing Eqs. (17, 18) we conclude that the potential 
Eq. (12) overestimates the dephasing effect and gives too 
large a value of a for N^Q^ <C N. Close to the cut-off 
p(X) oc VA — a, A i=a a independently of the exact form of 
W{\) and we conclude that for p(T) w — T,T w 

T + effects of dephasing and asymmetric contacts are sim- 
ilar. 

Even though p(X) does not behave as the coherent 
distribution po(A) = N/(2ttV~X), it should converge to 
po(A) as a — > 0. While its exact form depends on 
W(X), we can factorize it using the unknown functions 
f,g, which are smooth on a small scale <C 1, as p(X) = 
(7V/27r)VA _ ^/(a/A). 9 (A). Keeping A/a > 1 fixed, from 
the asymptotic po(A) we find g(X) = 1/A and /(0) = 1. 
For a 7^ the distributions equalize, p(X) = po(X), when 
yjl — a/Xf(a/X) = 1. Obviously, the second solution to 
this equation should exist at a/A < 1 (for example, for 
p(A) given by Eq. (16) the root is a/A = V3/2) and as 
a — > 0, this solution moves, A — > 0. In T-variables not a 
divergence po(7~) oc l/vl — 7~, but rather a peak of p(T) 
is then expected, p(T) — l/y^/Erh, 1 - 7" ~ j<j,/E Th 
and the bulk of p(T) is only slightly perturbed compared 
to po{T). Schematically this distribution is presented in 
the inset to Fig. M2. 

The important result (18) shows that the dephasing- 
induced gap for transmissions close to T = 1 survives the 
limit N » 1, 1 — T + w j^/2TrE Th = ra/r^. The shape 
of p(T) cx (T+ - T)^ 2 ,T « T+ for AT > 1 should be 
contrasted with p(T) cx exp[-/3iV ^/2(l - T)],T w 1 
for TV = 2. We conclude that lowering temperature T, 
7^ diminishes and thus opens conducting channels. In a 
normal system the widening of p{T) and the growth in 
VarC/ usually come together. In some models, like the 



Fabry-Perrot interferometer with T± given in Eq. (7) , or 
a weakly dephased quantum dot with 7± presented in Eq. 
(18), Var^ is defined by (T+ - 71 ) 2 only, see Eq. (14). 
Even though the connection of Var Q to T± is generally 
much more complicated [18], we suggest that a wider dis- 
tribution is a signature of enhanced Y&vQ. Since the 
normal current is only proportional to T, I{T) cx VT , 
this enhancement is hardly noticeable in a dephased dot 
at small j4,/E T h < 1, see Eq. (14). However, the SNS 
current is exponentially dominated by perfect channels, 
I(T) cx cxp(— 7rA(l — T)/eV), and a weak dephasing 
j^/ETh *C 1 remains important for SNS transport fluc- 
tuations. In a realistic model of our SNS experiment 
both dephasing and imperfect contacts should be taken 
into account. 
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